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1 Introduction 



This paper is concerned with the elastoplastic behavior of and moisture diffusion through 
a nanocomposite consisting of a vinyl ester resin matrix filled with montmorillonite clay 
particles. 

The choice of vinyl ester resin is explained by the fact that this thermosetting polymer 
is widely used as a matrix for glass-reinforced polymer composites employed for construc- 
tion and repair of bridges and other civil structures (Valea et al., 1998; Patel et al., 
1999). 

Montmorillonite (MMT) is a commonly used clay with a layered structure that is 
constructed of two tetrahedral sheets of silica surrounding an octahedral sheet of alumina 
or magnesia (Masenelli-Varlot et al., 2002). The layers (with thickness of 1 nm) are stacked 
by weak dipole forces, while the galleries between the layers are occupied by metal cations. 
Processing of clay particles for preparation of an intercalated nanocomposite consists of 
two stages: 

1. The growth of the basal (interlayer) spacing by exchanging the metal cations with an 
intercalating reagent (conventionally, alkylammonium ions) . This process facilitates 
intercalation of polymer chains into the galleries between the clay layers (Chang et 
al., 2002). 

2. Treatment of intercalated particles by a compatibilizer (a swelling agent, a functional 
oligomer, or a low molecular weight polymer) to improve miscibility between the 
organics- modified silicate layers and the host matrix (Tyan et al., 2000; Ma et al., 
2001; Chang et al., 2002). 

Preparation and analysis of micro-structure and physical properties of nanocomposites 
with polymeric matrices filled with MMT clay have been a focus of attention in the past 
couple of years. It was found that polymer-clay nanocomposites possess many desirable 
properties, such as increased moduli and strength, reduced gas permeabilities, decreased 
adsorption of organic liquids, etc. (Kim et al., 2002). However, no detailed analysis of 
interaction between molecules of a host matrix and the silica sheets has been performed 
so far. 

The objective of the present study is to investigate how the presence of intercalated 
MMT clay affects the mobility of polymeric chains. We focus on two main processes, 
where this mobility is revealed: 

1. elastoplasticity (which is attributed to sliding of junctions in the polymeric network 
with respect to their reference positions in the bulk material), 

2. moisture diffusion (whose rate is governed by the ability of penetrant molecules to 
move through the host network). 

In civil engineering applications, these two processes are of essential importance for the 
prediction of durability and performance of constructions. 

To evaluate the effect of nano-particles on the mobility of macromolecules, we focus 
on a model nanocomposite where no compatibilizer is used. Because of low miscibility 
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between the clay particles and the vinyl ester resin, we do not expect that the presence of 
intercalated clay particles noticeably improves the mechanical and transport properties 
of a nanocomposite compared to those of the neat resin (Ma et al., 2001). Our aim 
is to demonstrate that the concentration of nano-particles affects the parameters that 
characterize mobility of polymeric chains in mechanical and diffusion tests in similar 
ways. 

In the past two decades, a number of constitutive models have been proposed for the 
elastoplastic and viscoplastic behavior of glassy polymers, see, e.g., Boyce et al. (1988), 
Mangion et al. (1992), Hasan and Boyce (1995), David et al. (1997), Drozdov (2001) and 
the references therein. Three shortcomings of these models should be mentioned: 

• the stress-strain relations were developed to correctly describe the mechanical be- 
havior of polymers near the yield point (which is not the case for vinyl ester resin 
where fracture occurs earlier than yielding, see Figures 1 to 5 below), 

• the constitutive equations contain a number of adjustable parameters which cannot 
be determined with a high level of accuracy by fitting data in uniaxial tensile tests, 

• the material constants are not explicitly expressed in terms of the physical quantities 
that describe molecular mobility. 

To assess constrains imposed by the clay particles on the mobility of macromolecules in 
a host matrix, we develop a simple model that contains only three adjustable parameters 
with a transparent physical meaning. Only one of these parameters reflects mobility of 
chains, and it can be easily found by fitting stress-strain curves in uniaxial tensile tests 
with small strains. 

Diffusion of moisture in polymer composites has attracted substantial attention in the 
past decade, see, e.g., Cai and Weitsman (1994), Bavisi et al. (1996), Vanlandingham et 
al. (1999), Marcovich et al. (1999), Roy et al. (2000), Chen et al. (2001), Roy et al. 
(2001), Uschitsky and Suhir (2001), to mention a few. It was demonstrated that diffusion 
of penetrants in several glassy polymers is anomalous (non-Fickian), with the mass gain 
proportional to some power of time that differs from 0.5 (the latter is characteristic for 
the Fickian diffusion). 

Two approaches were recently proposed to model the anomalous diffusion. According 
to the first (the fractional dynamics concept), conventional derivatives with respect to 
time and spatial variables in the diffusion equation are replaced by fractional derivatives, 
see Metzler and Klafter (2000) and the references therein. According to the other (Roy at 
al., 2000), a constant coefficient of diffusion is replaced by a decreasing function of time 
(by analogy with a relaxation modulus for a viscoelastic solid). 

Our experimental data (see Figures 6 to 10 below) reveal that moisture diffusion in 
the neat vinyl ester resin is practically Fickian, while the anomalous transport of water 
molecules is observed with an increase in the clay concentration. This phenomenon may 
be explained by the fact that the water molecules are bounded to hydrophilic surfaces of 
clay layers, where they become immobilized. A two-stage diffusion process with immobi- 
lization of penetrant molecules has been previously studied by Carter and Kibler (1978) 
and Gurtin and Yatomi (1979). Unlike these works, we assume that bounded molecules 
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cannot leave the sites where they were immobihzed. This allows the number of adjustable 
parameters to be substantially reduced. The diffusion process is governed by differential 
equations with only three material constants that are found by matching the water uptake 
curves plotted in Figures 6 to 10. 

The exposition is organized as follows. Observations in uniaxial tensile tests and in 
moisture diffusion tests are reported in Section 2. Constitutive equations for the elasto- 
plastic response are derived in Section 3 by using the laws of thermodynamics. Adjustable 
parameters in these relations are determined in Section 4. A model for moisture diffusion 
is developed in Section 5 and its material constants are found in Section 6. A brief dis- 
cussion of our findings is presented in Section 7. Some concluding remarks are formulated 
in Section 8. 

2 Experimental procedure 

2.1 Preparation of samples 

The polymer used was DERAKANE 411-350 epoxy vinyl ester resin (Dow Chemical 
Co.) containing 45 wt.-% of dissolved styrene. To cure the resin at room temperature, 
it was mixed with 0.5 wt.-% of 6% cobalt naphthenate catalyst (Sigma Aldrich Co.). 
Additionally, 0.05% of 99% N,N dimethyl aniline (Lancaster Synthesis, Pelham, NH, 
USA) was used as an accelerator, and 1.5% of methyl ethyl ketone peroxide with 9% of 
active oxygen (Sigma Aldrich Co.) as an initiator. 

A montmorillonite clay Cloisite lOA (Southern Clay Products Inc., Gonzales, TX, 
USA) treated with benzyl (hydrogenated tallow alkyl) dimethyl quaternary ammonium 
chloride was used as received. Neat resin coupons were cast by pouring the reaction 
mixture into Teflon molds with the dimensions 50 x 12.5 mm and the thickness ranging 
from 0.2 to 0.6 mm. The organically treated clay was added to the liquid resin and 
manually stirred. The mixture was then degassed in a vacuum oven to remove air bubbles. 
Afterwards, the catalyst, initiator and accelerator were added. The mixture was allowed 
to cure at room temperature for 24 hours, and it was subsequently post cured in an oven 
for 3 hours at 90 °C. 

TEM (transmission electron microscopy) micrographs of the nanocomposite samples 
show that the clay is reasonably uniformly distributed, and the silica sheets are randomly 
oriented (Shah et al., 2002). The clay exists as expanded aggregates made up of 2 to 
10 sheets with the interlayer distance ranging from 4.4 to 5.0 nm, which confirms the 
formation of an intercalated nanocomposite. 

2.2 Mechanical tests 

Dumbbell specimens for mechanical tests were cut following ASTM D638 specification. 

Tensile tests were performed at ambient temperature by using a 100 kN Instron machine 
model 8501 with a cross-head speed of 0.254 mm/min. This strain rate ensures prac- 
tically isothermal loading conditions. The strain was measured independently using a 
strain gauge affixed to the mid-point of the specimen. The tensile force was measured by 
a standard load cell. The engineering stress, cr, was determined as the ratio of the tensile 
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force to the cross-sectional area of a specimen in the stress-free state. For any concen- 
tration of filler in a nanocomposite, four specimens were tested. Typical stress-strain 
diagrams are presented in Figures 1 to 5. These figures demonstrate that the response 
of the neat resin is strongly nonlinear, but the nonlinearity of the stress-strain curves is 
reduced with an increase in the filler content. 

DSC (differential scanning calorimetry) measurements show that the glass transition 
temperature of the neat vinyl ester resin is about 97 °C, and it increases with concentration 
of clay reaching 117 °C at = 5.0 wt.-% (Shah et al., 2002). Because the material at room 
temperature is far below its glass transition point, we disregard its viscoclastic response 
and attribute the nonlinearity of the stress-strain curves to the elastoplastic behavior. 



2.3 Diffusion tests 

Diffusion tests were performed by immersing samples with rectangular cross-section in 
distilled water at room temperature (T = 25 °C). The samples had a dry mass ranging 
from 120 to 400 mg. They were stored in a controlled humidity chamber, and contained 
0.05 ± 0.005 wt.-% water at the beginning of diffusion tests. The samples were periodically 
removed, blotted dry with a lint-free tissue, weighed and re-immersed in the water. A 
typical experiment lasted ten days, and, on the first day, readings were taken as frequently 
as every 15 min. The balance used had an accuracy of 1 mg, and 3 to 5 replicate runs 
were carried out for a given set of condition. Typical dependences of the relative moisture 
uptake, $, versus the reduced time, t, are plotted in Figures 6 to 10. The reduced time, 
t, is given by 

t 

where t is time elapsed from the beginning of a test, and 21 is thickness of a sample. The 
relative moisture uptake, $, is defined as 



*■= 775. (1) 



K[oo) 

where R(t) is the mass gain at the reduced time t and R{oo) is the maximal mass gain. 
Observations on samples with various thicknesses (not presented) show that the mass gain 
curves plotted in the reduced coordinates ($, t) are superposed fairly well. 

Figures 6 to 10 demonstrate that moisture diffusion in the neat resin is roughly Fickian, 
but it becomes strongly anomalous with an increase in the filler content. 

The aim of the remaining part of this work is to develop constitutive equations for the 
elastoplastic behavior and moisture diffusion in a nanocomposite and to find adjustable 
parameters in these relations by matching the experimental data plotted in Figures 1 to 
10. 



3 A model for the elastoplastic response 

To simplify the analysis, we model a nanocomposite as an equivalent network of macro- 
molecules bridged by junctions (entanglements, physical cross-links and particles of filler). 
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Because the viscoelastic behavior of the nano composite is disregarded, the network is 
treated as permanent (we suppose that active chains do not separate from junctions and 
danghng chains do not merge with the netwotk within the experimental time-scale). To 
describe the elastoplastic behavior of the network, we assume that junctions slide with 
respect to their reference positions in the bulk material. For uniaxial deformation, sliding 
of junctions is determined by a plastic strain ep. Adopting the conventional hypothe- 
sis that the macro-strain, e, is transmitted to all chains in the network by surrounding 
macromolecules, we write 

e = ee ep, (3) 

where is an elastic strain. We accept the mean-field approach and treat and ep as 
average elastic and plastic strains per chain. 

It is assumed that under active loading, the rate of changes in the strain, ep, with 
time, t, is proportional to the rate of changes in the macro-strain e, 

(4) 

where the coefficient of proportionality, is a function of the clastic strain e^- 
The following restrictions are imposed on the function if^e^): 

1. The coefficient of proportionality in Eq. (4) vanishes at the zero elastic strain, 
(^(0) = 0. 

2. The function </? monotonically increases with ee and tends to some constant a e (0, 1) 
for relatively large elastic strains, 

lim (fiee) = a. 

The constant a determines the rate of shding in junctions at the stage of a developed 
plastic flow. 

To approximate experimental data, we use the function 

(^(ee) = a [l - exp (-—)], (5) 



which satisfies the above conditions. This function is determined by two adjustable pa- 
rameters, a and e, where a is the hmiting value of </? at "large" elastic strains, and e > 
is a constant that characterizes how "large" an elastic strain, ee, is. 

A chain is modelled as a linear elastic solid with the mechanical energy 

1 2 

where is an average rigidity per chain. Multiplying the energy, w, by the number of 
chains per unit volume, no, we find the strain energy density 

W = \Eel (6) 
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where E = uriQ is an elastic modulus. It is worth noting that the elastic modulus, may 
differ from the Young modulus. The latter is conventionally defined as the tangent of 
the angle between the tangent straightline to the stress-strain curve at small strains and 
the horizontal axis. These two quantities coincide when the elastic strain, ee, is equal to 
the macro-strain, e. For the elasto-plastic response of a nanocomposite, when the stress- 
strain diagram substantially difi^ers from a straightline, E exceeds the Young modulus 
determined by the traditional method. 

For isothermal uniaxial deformation, the Clausius-Duhem inequality reads, 

^ , , dW , , , .de , . 

= -_(,) + .(*)-(*) >0, 

where Q is internal dissipation per unit volume. Substition of Eqs. (3), (4) and (6) into 
this formula results in 



Q{t)^[a{t)-Eeeit){l-if{e,{t))) 



de , . 



Assuming the expression in square brackets to vanish, we arrive at the stress-strain rela- 
tion 



(7) 



a(t) = Ee,(t) 

Constitutive equations (3) to (5) and (7) are determined by 3 adjustable parameters: 

1. the elastic modulus E, 

2. the rate of developed plastic fiow a, 

3. the strain, e, that characterizes transition to the steady-state plastic deformation. 
To find these quantities, we match experimental data depicted in Figures 1 to 5. 



4 Fitting observations in mechanical tests 

It follows from Eqs. (3), (4), (5) and (7) that in a uniaxial tensile test, the stress is given 
by 

aie) = ^(e - ep){l - a[l - exp(-i-^)] }, (8) 
where the plastic strain, ep, satisfies the nonlinear differential equation 



dep 
de 



exp 



(9) 



with the initial condition ep(0) = 0. 

We begin with matching the experimental data in a test on the neat vinyl ester resin, 
see Figure 1. To find the material constants, E, a and e, we fix some intervals [0, flmax] and 
[0,emax], where the "best-fit" parameters a and e are assumed to be located, and divide 
these intervals into / subintervals by the points = iAa and ej = jAe = !,...,!) 
with Aa — a-aiax/I and Ae = emax/I- For any pair, {ai,ej}, we integrate Eq. (9) 
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numerically (with the step Ae = 1.0 • 10~^) by the Runge-Kutta method. Given {a^, cj}, 
the elastic modulus E — E{i, j) is found by the least-squares technique from the condition 
of minimum of the function 




where the sum is calculated over all experimental points, e^, depicted in Figure 1, CTexp 
is the stress measured in a tensile test, and (Xnum is given by Eq. (8). The "best-fit" 
parameters a and e minimize the function F on the set {a,, ej {i, j — 1, . . . , I)}. Fitting 
the observations results in e = 2.11 ■ 10~^. 

Wc fix this value of e and proceed with matching observations in tensile tests on other 
specimens by using only two material constants, E and a, which arc determined by the 
above algorithm. The average values of E and a (over 4 samples) are depicted in Figures 
11 and 12, where the vertical bars stand for the standard deviations of these parameters. 

Afterwards, we match experimental data in tensile tests on specimens with various 
contents of filler, u. The same algorithm of fitting is employed with the fixed value 
e = 2.11-10"^. Figures 1 to 5 demonstrate fair agreement between the observations and the 
results of numerical simulation. The elastic modulus, E, and the rate of developed plastic 
flow, a, are plotted versus the filler content, in Figures 11 and 12. The experimental 
data are approximated by the phenomenological equations 



where the coefficions Ej and aj {j — 0, 1) are found by the least-squares technique. 

5 A model for the moisture diffusion 

A sample is treated as a rectilinear plate with thickness 21. We introduce Cartesian 
coordinates {x, y, z}, where the axis x is perpendicular to the middle plane of the plate, 
and the axes y and z lie in the middle plane. Length and width of the plate are assumed to 
substantially exceed its thickness, which implies that the moisture concentration depends 
on the coordinate x only. 

After the plate is immersed into water, three processes occur in the composite material: 

1. sorption of water molecules to the plate faces from the surrounding, 

2. diffusion of the penetrant into the plate, 

3. adsorption of water molecules on the surfaces of clay layers, where these molecules 
become immobilized. 

It is conventionally accepted that the rate of sorption in glassy polymers noticeably ex- 
ceeds the rate of diffusion (Chen et al., 2001), which implies that the sorption equilibrium 
is rapidly established. The equilibrium condition is given by 



E 



Eq — EiU, 



a = Qq — ttiu, 



(10) 



n{t, x) 



x=±l 



— n 



o 



(11) 
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where n is the moisture concentration at time t at point x [the number of water molecules 
in the polymeric matrix occupying a volume with the unit area in the plane (y, z) and 
the thickness dx reads n{t, x)dx ], and n° is the equilibrium moisture concentration in the 
matrix on the faces of the plate. 

Diffusion of the penetrant through the matrix is described by the mass conservation 

law 

dn _dJ_ _dni 

dt~dxdt' ^ ^ 

where J{t,x) stands for the mass flux, and ni{t,x) denotes the concentration of water 
molecules immibilized at the surfaces of clay layers. 
The mass flux obeys the Fick equation 

. = (13) 

where D is diffusivity. Equation (13) is tantamount to the assumption that diffusion 
occurs through the polymeric matrix only. However, the coefficient of diffusion, D, is 
assumed to depend on the volume content of filler, because the presence of particles 
results in a decay in the molecular mobility of polymeric molecules which provides the 
driving force for diffusion of penetrants. 

Adsorption of water molecules on the surfaces of filler is determined by the first-order 
kinetic equation 

^^kn{nl-n,). (14) 

According to Eq. (14), the rate of adsorption is proportional to the concentration of 
the penetrant in the matrix, n, and to the current number of "unoccupied sites" on the 
surfaces of filler, — ni, where nl is the total number of sites where water molecules can 
be immobilized. In the general case, the rate of adsorption, k, and the maximal concen- 
tration of unoccupied sites, nl, are functions of the clay content, because the presence of 
nanoparticles affect mobility of chains in the polymeric matrix, and, as a consequence, 
their chemical potential. 

Neglecting the moisture content in a sample before testing, we adopt the following 
initial conditions: 

nit,x)l^=0, mit,x)l^=0. (15) 

Equations (12) to (14) with initial conditions (15) and boundary condition (11) uniquelly 
determine the diffusion process. The moisture mass gain (per unit mass of the sample) is 
given by 

R{t) = —— I n{t,x) + ni{t,x) dx, (16) 



2lp 

where n is the average mass of a penetrant molecule, and p is mass density of the matrix. 

Our aim now is to transform Eqs. (11) to (16). For this purpose, we introduce the 
dimensionless coordinate, x — x/l, and the reduced moisture concentrations, c = n/n° 
and Ci = ni/n°. Combining Eqs. (1), (12) and (13), we arrive at the mass flux equation 

dc f^d'^c dc 
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where 

In the new notation, Eq. (14) reads 

dci 



^^Kc(C-c,), (19) 



where 

K = Akfn°, C = ^. (20) 



Initial conditions (15) remain unchanged, 

c(0,x) = 0, ci(0,x) = 0. (21) 

Bearing in mind that the functions c{t, x) and ci{t, x) are even functions of x, we present 
boundary condition (11) in the form 

Oc - - 

— (f,0) = 0, c(i,l) = l. (22) 

Finally, Eq. (16) is given by 

R(t) = -n° r \c(i, x) + ci it, x)] dx. (23) 

p Jo L J 

With an increase in time, t, the solutions of Eqs. (17) and (19) tend to the steady-state 
solutions 

c(oo,a;) = 1, Ci(oo,a;) = C. (24) 
It follows from Eqs. (23) and (24) that the maximal moisture uptake reads 

i?(oo) = ^n°(l + C), (25) 
whereas the relative mass gain, see Eq. (2), is given by 

1 + G JO L 

Equations (17), (19) and (26) with initial conditions (21) and boundary conditions (22) 
are determined by 3 adjustable parameters: 

1. the reduced diffusivity D, 

2. the rate of moisture adsorption on nanoparticles 

3. the reduced concentration of unoccupied sites on the surfaces of filler C. 
These constants are found by fitting the experimental data in diffusion tests. 
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6 Fitting observations in diffusion tests 



We begin with matching the experimental data for the neat vinyl resin depicted in Figure 
6. Because in the absence of filler, moisture is not adsorbed on the surfaces of particles, 
we set C = in Eq. (17). Equations (17), (21) and (22) are reduced to the conventional 
boundary problem for the diffusion equation 

= ^(^"1) = !' ^(0'^) = 0' (27) 

which is resolved by the finite- difference method with an explicit algorithm. Given a 
coefficient of diffusion, Z), we divide the interval [0, 1] into M subintervals by points 
Xm = mAx {m = 0, 1, ... , M) with Ax = 1/M, introduce discrete time in = nAt, and 
replace Eq. (27) by its finite-difference approximation 



c{0,Xm) = (m = 0,1,..., M-1), c{0,xm) = 1, 

At_ 

Ax"' 

+c{in,Xm-i)\ (m = 1, . . . ,M - 1), 
c(t„+i,Xo) = c(t„+i,xi), c(t„+i,XM) = l (n = l,2,...). (28) 

To calculate the relative mass gain, the integral in Eq. (26) is approximated by the Euler 
formula 

M-1 

$(?„) = Ax ^ cit ). (29) 

m=0 

The coefficient D is determined by the following procedure. We fix some interval 
[-Dmin, -Dmax], where the "best-fit" value of D is supposed to be located, and divide this 
interval into I sub- intervals by the points Di — Dmin + iAD [i — 0,1,...,/), where 
AD = (Z^max — Dnnn)/I- For any Di, Eq. (28) is solved numerically (with the steps 
Ax = 0.05 and At = 0.005 that ensure the stability of the numerical algorithm), and the 
function $(tn) is given by Eq. (29). The "best-fit" value of D is found from the condition 
of minimum of the function 



F = ^[*exp(tfe) -*num(^fe) 



where the sum is calculated over all point ik depicted in Figure 6, $cxp(^fe) is the relative 
mass gain measured in the test, and $num(^^feO is determined from Eq. (29). After the "best- 
fit" diffusivity, Di, is found, the procedure is repeated for the new interval [-Dj-i, Di+i] to 
guarantee an acceptable accuracy of matching observations. Figure 6 demonstrates fair 
agreement between the experimental data and the results of numerical simulation. 

To ensure that the numerical results are independent of the choice of the steps of 
integration. Ax and At, we repeat integration of Eqs. (28) and (29) with the "best-fit" 
parameter D and the new values Ax = 0.01 and At = 0.0005. No difference was observed 
between the curve plotted in Figure 6 and the curve obtained by the numerical integration 
of Eqs. (28) and (29) with the decreased steps in time and the spatial coordinate. 
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The above procedure is repeated to match observations in the tests on 4 different 
samples. The average value of diffusivity is presented in Figure 13, where the vertical 
bars stand for the standard deviation. 

It is worth noting that the algorithm of fitting experimental data used in this work 
differs from the conventional procedure employed by Shah et al. (2002), according to 
which the coefficient of diffusion was determined by matching the initial parts of the mass 
gain diagrams (where the function $(t) is approximately linear). 

The same numerical algorithm is employed to match experimental data on polymeric 
nanocomposites with various amounts of filler. The difference in the treatment of ob- 
servations consists in the following: (i) Eq. (28) arc replaced by the finite-difference 
approximation of Eqs. (17) and (19); (ii) the "best-fit" values are simultaneously deter- 
mined for 3 adjustable parameters: D, C and K. Figures 7 to 10 show good agreement 
between the experimental data and the results of numerical analysis. The average values 
of D [this quantity is found from Eq. (18)], C and K (over 4 tests) are plotted versus 
the filler fraction, v, in Figures 13 to 15. The dependences of these parameters on the 
concentration of filler are approximated by the linear equations 

D^Do-Diu, C^Ciu, K^Ko-Kiu, (30) 

where the coefficions Dj, Cj and Kj (j = 0, 1) are found by the least-squares technique. 

It follows from Eq. (25) that the maximal moisture uptake per unit mass of the matrix, 
-Rpoi = Kn°/p, is given by 

_ Rjoo) . . 

- YTc' ^ ' 

Using experimental data for i?(oo) and the values of C determined by matching observa- 
tions, we calculate Rpoi from Eq. (31) and plot this parameter versus the clay content in 
Figure 16. 



7 Discussion 



Figure 11 demonstrates that the elastic modulus, E, decreases with the clay content, u, 
in the region of small concentrations of filler {u < 1.0 wt.-%), and remains practically 
constant when u exceeds 1.0 wt.-%. According to Figure 16, the filler fraction, u, affects 
the maximal moisture uptake, i?poi, in a similar way. These dependences of the parameters 
E and i?poi seem quite natural, because these quantities characterize the macro-response 
of a nanocomposite in tensile and diffusion tests. They may be explained by a progressive 
decrease in the degree of cross-linking with the growth of clay concentration (Bharadwaj 
et al., 2002), which results in a decrease in the number of active chains (whose ends are 
linked to nearby junctions) compared to the neat vinyl ester resin. This implies that the 
elastic modulus (proportional to the number of active chains) is reduced and the maximal 
weight gain (proportional to the free volume per unit active chains) is increased. Figures 
11 and 16 show that the concentration of chains reaches its limiting value rather rapidly 
(when the filler content is about Uc = 1 wt.-%), and an increase in the clay concentration 
above this critical value does not affect the cross-linking process. 
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Figure 12 demonstrates that the rate of developed plastic flow, a, strongly decreases 
with 1/ in the region below the critical concentration of filler, v < Vc, and it weakly 
decreases with u when the clay content exceeds its critical value. Figures 13 and 15 reveal 
that the coefficient of diffusion, D, and the rate of adsorption of water molecules on 
surfaces of nano-layers, demonstrate similar dependences on the filler fraction. This 
behavior of the quantities a, D and K appears to be in agreement with our physical 
intuition, because all these parameters reflect mobility of polymeric chains. According to 
Figures 12, 13 and 15, this mobility is drastically reduced with v in the sub-critical region, 
and it proceeds to decrease (however, rather weakly) when the flUer content exceeds the 
critical concentration v^- The conclusion that the mobility of chains in the polymeric 
matrix decreases with the clay content is also conflrmed by DSC measurements (Shah et 
al., 2002) which reveal the the glass transition temperature monotonically increases with 
I/. 

Figure 14 reveals that the parameter C linearly increases with v. It follows from Figure 
16 and Eqs. (25) and (31) that the maximal moisture content in the polymeric matrix, 
n°, is independent of the clay fraction. This result together with Eq. (20) and Figure 14 
implies that number of sites on the clay layers where water molecules are immobilized, n^, 
linearly grows with the flUer concentration. The latter assertion is in agreement with the 
main assumptions of our model (the total number of sites where the penetrant molecules 
are immobilized equals the product of the average number of these sites per particle by the 
number of clay particles) , and it may be treated as a conflrmation of the kinetic equations 
for the anomalous moisture diffusion. 

The similarity in this behavior of the quantities a, D and K as functions of the clay 
concentration, provides a way to substantially reduce the number of (time-consuming) 
moisture diffusion experiments. If the parameters a and D decrease with v in the same 
fashion (our analysis of observations on vinyl ester resin nanocomposite with MMT par- 
ticles treated with vinyl benzine trimethyl ammonium chloride conflrms this hypothesis), 
the function a{u) can be determined with a high level of accuracy in standard tensile 
tests, whereas only a few moisture diffusion tests are necessary to flnd the dependence of 
diffusivity on the clay content, D — D{v). 

8 Concluding remarks 

Mechanical and diffusion tests have been performed on a nanocomposite with vinyl ester 
resin matrix and montmorillonite clay flUer at room temperature. Experimental data in 
tensile tests demonstrate that the mechanical behavior of the neat vinyl ester resin is 
strongly elastoplastic, whereas an increase in the clay content results in a decrease in the 
plastic strain (due to the constrains on mobility of junctions between polymeric chains 
imposed by nano-particles) . Observations in diffusion tests show that moisture diffusion 
in the neat resin is Fickian, whereas it becomes noticeably anomalous (non-Fickian) with 
the growth of the clay content. This transition is ascribed by adsorbtion (immobilization) 
of the penetrant molecules on surfaces of hydrophilic clay layers. 

Constitutive models have been proposed to describe the elastoplastic response and 
anomalous moisture diffusion in a nanocomposite. Both models are determined by 3 
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adjustable parameters that are found by fitting the experimental data. Fair agreement is 
demonstrated between the observations and the results of numerical simulation. 
The following conclusions are drawn: 



1. A critical concentration of clay is found, i^c ~ 1 wt.-%, corresponding to both 
mechanical and diffusion tests. The filler content, v, strongly affects the material 
constants in the sub-critical region, and its influence becomes rather weak in the 
post-critical domain. 

2. The quantities independent of molecular mobility, such as the elastic modulus, 
and the maximal moisture uptake by the polymeric matrix, R^oh are not affected by 
the concentration of clay particles when the filler content exceeds its critical value. 

3. The quantities reflecting mobility of polymeric chains, such as the rate of developed 
plastic flow, a, the coefficient of diffusion, and the rate of moisture adsorption on 
nano-particles, A;, strongly decrease with the filler content in the sub-critical region, 
and proceed to decrease (with smaller rates) ai u > Uf.. 
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Figure 1: The stress a MPa versus strain e for neat vinyl ester resin. Circles: experimental 
data. Solid line: results of numerical simulation 



17 



40.0 



a 



0.0 




0.0 e 0.01 



Figure 2: The stress a MPa versus strain e for nanocomposite with v — 0.5 wt.-%. Circles: 
experimental data. Solid line: results of numerical simulation 
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Figure 3: The stress a MPa versus strain e for nanocomposite with v — 1.0 wt.-%. Circles: 
experimental data. Solid line: results of numerical simulation 
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Figure 4: The stress a MPa versus strain e for nanocomposite with v — 2.5 wt.-%. Circles: 
experimental data. Solid line: results of numerical simulation 
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Figure 5: The stress a MPa versus strain e for nanocomposite with v — 5.0 wt.-%. Circles: 
experimental data. Solid line: results of numerical simulation 
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Figure 6: The relative moisture uptake $ versus the reduced time i (hz/mm) for neat 
vinyl ester resin. Circles: experimental data (21 = 0.56726 mm). Solid line: results of 
numerical simulation 
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Figure 7: The relative moisture uptake $ versus the reduced time i (h^/mm) for nanocom- 
posite with u = 0.5 wt.-%. Circles: experimental data {21 — 0.25373 mm). Solid line: 
results of numerical simulation 
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Figure 8: The relative moisture uptake $ versus the reduced time i (ha /mm) for nanocom- 
posite with u = 1.0 wt.-%. Circles: experimental data {21 — 0.20108 mm). Solid line: 
results of numerical simulation 
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Figure 9: The relative moisture uptake $ versus the reduced time t (h^ /mm) for nanocom- 
posite with u = 2.5 wt.-%. Circles: experimental data {21 — 0.22013 mm). Solid line: 
results of numerical simulation 
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Figure 10: The relative moisture uptake $ versus the reduced time i (h^/mm) for 
nanocomposite with u = 5.0 wt.-%. Circles: experimental data {21 — 0.2032 mm). Solid 
line: results of numerical simulation 
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Figure 11: The elastic modulus E GPa versus the clay concentration u wt.-%. Circles: 
treatment of observations. Solid lines: approximation of the experimental data by Eq. 
(10). Curve 1: Eq = 4.83, Ei = 0.85; curve 2: Eq = 3.97, Ei = 0.0 
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Figure 12: The rate of developed plastic flow a versus the clay concentration v wt.-%. 
Circles: treatment of observations. Solid lines: approximation of the experimental data 
by Eq. (10). Curve 1: Oq = 0.22, ai = 0.19; curve 2: oq = 0.04, oi = 0.01 
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Figure 13: Diffusivity D ■ 10 ^ mm^/s versus the clay concentration u wt.-%. Circles: 
treatment of observations. Solid lines: approximation of the experimental data by Eq. 
(31). Curve 1: Dq = 9.01, Di = 3.91; curve 2: Dq = 6.11, Di = 0.24 
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Figure 14: The dimensionless parameter C versus the clay concentration v wt.-%. Circles: 
treatment of observations. Solid line: approximation of the experimental data by Eq. (31) 
with Ci = 0.19 
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Figure 15: The adsorption rate iT- 10 ^ versus the clay concentration z/ wt.-%. Circles: 
treatment of observations. Solid lines: approximation of the experimental data by Eq. 
(31). Curve 1: Kq = 14.68, = 12.85; curve 2: Kq = 2.14, = 0.19 
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Figure 16: The maximal moisture uptake per unit mass of the matrix i?poi % versus the 
clay concentration i> wt.-%. Circles: treatment of observations. Solid line: approximation 
of the experimental data by the constant Rpoi = 0.75 % 
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